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CHAPTER  I 


INTRODUCTION 

Logistical  support  of  an  operational  weapon  system 
requires  that  spares  be  available  to  replace  components 
that  fail.  The  cost  of  the  complex  equipment  used  in 
modern  weapon  systems  dictates  that  the  number  of  spares 
in  the  logistics  system  be  kept  to  a minimum.  Determina- 
tion of  minimum  spares  requirements  depends,  upon  an  accur- 
ate prediction  of  expected  failures. 

Background 

' 

The  Air  Force  distinguishes  between  initial  and 
replenishment  spares.  Usually,  initial  spares  support 
a weapon  system  from  the  time  of  its  preliminary  opera- 
tional capability  through  the  item  lead  time  plus  throe 
months.  In  no  case,  is  the  initial  support  period  less 
than  twelve  months  <2:p.l-l).  Replenishment  spares  sup- 
port the  weapon  system  thereafter. 

Spares  are  classified  as  either  recoverable  or 
expense  items.  Recoverable  items  are  repaired  and 
returned  to  service  if  they  fail;  expense  items  are  thrown 
away.  Expense  item  replenishment  requirements  are  com- 
puted under  Economic  Order  Quantity  procedures  (5) . More 
complex  procedures  for  computing  recoverable  item 
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replenishment  requirements  are  embodied  in  the  Recoverable 
Consumption  Item  Requirements  System,  designated  D041 
(3 : p.1-1) . 

The  D041  System.  The  D041  system  operates  at 
each  of  the  Air  Force  Logistics  Command's  five  Air  Logis- 
tics Centers.  A segment  of  the  system  also  operates 
at  HQ  AFLC.  D041  is  designed  to  accomplish  tho  following: 

a.  Compute  requirements  for  recoverable  . . . 
items  .... 

d.  Perform  tho  routine  clerical,  mathematical, 
and  statistical  workload  involved  in  computing  recov- 
erable item  requirements. 

c.  Forecast  gross  and  not  requirements  using  past 
and  future  programs,  usage  history,  and  asset  informa- 
tion maintained  within  this  system. 

d.  Produce  reports  for  management  evaluation  and 
action • 

e.  Produce  information  for  other  automated  data 
systems  [ 3 :p . 2—1] . 

Recoverable  item  requirements  are  computed  by  the 
D041  system  once  every  quarter  of  tho  fiscal  year.  The 
first  computation  of  tho  fiscal  year  is  used  to  develop 
apportionment  roquosts  and  tho  third,  to  develop  budget 
requests.  All  four  computations  are  used  to  identify 
items  requiring  logistics  actions,  such  as  buy,  repair, 
disposal,  or  termination  of  oxisting  procurement  actions 
(3:pp.l-l  to  1-2) . 

The  D041  system  incorporates  a complex  requirements 
computation  algorithm.  The  algorithm  first  multiplies 
a forecast  of  future  program  activity  by  a series  of  factors 
to  generate  a gross  requirement.  The  gross  requirement 
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consists  of  an  operating  requirement,  a base  stock  level, 
a depot  stock  level,  plus  several  miscellaneous  categories 
of  requirements  that  have  no  bearing  on  the  subject  of  this 
research.  The  next  step  is  to  apply  to  the  gross  i.yuire- 
ment  those  assets  forecast  to  be  available.  In  so  doing, 
a forecast  is  made  of  the  base  and  depot  level  repair  to 
be  accomplished  in  the  computation  period.  Requirements 
that  cannot  be  satisfied  by  the  inventory,  with  mainte- 
nance support,  must  be  supported  by  additional  procure- 
ment (3 :Ch.6)  . 

The  DO 41  computation  incorporates  a Variable 
Safety  Level  (VSL)  feature  based  on  the  METRIC  (Multi- 
Echelon  Technique  for  Recoverable  Item  Control)  model 
developed  for  the  Air  Force  by  The  RAND  Corporation.  A 
description  of  METRIC  can  be  found  in  a RAND  Memorandum, 
RM-5078-PR  (21) . Use  of  the  METRIC  model  requires  speci- 
fying the  probability  distribution  associated  wiv'-h  the 
demand  process  (17:7).  In  D041,  all  uemands  are  assumed  to 
occur  according  to  a Poisson  distribution  (15) . This 
research  examined  the  validity  of  the  assumption  as  it 
pertains  to  computing  requirements  for  specific  items. 

The  Poisson  Process.  A few  comments  on  the  Pois- 
son process  will  provide  the  statistical  concepts  under- 
lying the  problem  that  this  research  addresses.  "A  large 
class  of  situations  in  which  events  occur  randomly  is  the 
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Two  probability  distri- 


Poiaoon  process  [16:374-375]." 
butions,  the  Poisson  and  the  exponential,  are  associated 
with  a Poisson  process.  The  probabilities  for  the  number 
of  times  an  event  occurs  in  an  interval  of  time  follow  a 
Poisson  distribution;  while  the  probabilities  for  the  times 
between  occurrences  of  that  event  follow  an  exponential 
distribution  (16:375).  One  of  the  assumptions  that  must 
be  met  in  orde'  for  a Poisson  process  to  apply  is  that 
"the  process  rate  a must  remain  constant  for  the  entire 
duration  considered  [16:376]." 

Statement  of  the  Problem 

Recent  research  (0)  by  Captain?  Lowell  R.  Crowe 
and  Levi  D.  Lowman,  Jr.,  addressed  the  failure  patterns 
experienced  by  three  inertial  measurement  units  (IMUs) : 
the  FLIP  unit  used  on  the  C-5A  aircraft,  the  LN-15  unit 
used  on  the  D-52G/H  aircraft,  and  the  KT-73  unit  used  on 
the  n-7D/E  and  AC-130H  aircraft.  They  concluded  that  the 
boat  theoretical  distribution  for  describing  the  actual 
operating  time  to  failure  of  the  IMUs  was  not  the  expo- 
nential (8:44-48).  They  also  reported  indications  that 
the  IMU  failure  rates  may  have  changed  over  time  (8:71). 

Those  findings  call  into  question  whether  the 
failures  of  those  IMUs  can  be  legitimately  considered  to 
be  derived  from  a Poisson  process.  First,  of  course, 
the  applicability  of  the  exponential  distribution  itself 
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is  questioned,  although  the  evidence  does  not  conclusively 
eliminate  it  from  consideration.  The  inconclusiveness 
stems  from  the  decision  rule  used  to  determine  the  best 
distribution.  The  best  distribution  was  that  which 
described  the  most  samples  from  each  IMU  population.  In 
fact,  the  exponential  did  fit  some  of  the  samples 
(8:50-53).  Second,  the  assumption  of  a constant  process 
rate  was  undermined. 

Unless  the  failures  of  the  three  XMUs  can  be  con- 
sidered to  be  derived  from  a Poisson  process,  it  would  be 
inappropriate  to  use  the  Poisson  distribution  to  model 
demands  in  the  computation  of  requirements  for  these  items. 
The  problem  this  research  addressed  has  three  basic 
components.  First,  there  was  the  question  of  whether 
the  Poisson  process  does,  in  fact,  represent  the  failure 
pattern  for  the  three  IMUs.  Second,  t’.ero  was  the  ques- 
tion of  whether  substituting  the empii?cvl  failure  distri- 
bution for  the  IMUs  (or  another  appropriate  theoretical 
distribution)  in  the  computation  would  significantly 
change  the  computed  requirement.  Finally,  there  is  the 
possibility,  if  these  IMUs  do  not  exhibit  Poisson  charac- 
teristics, that  the  same  would  apply  to  other  IMUs. 

Approach  to  the  Research 

This  research  was  intended  to  pursue  the  first 
two  components  of  the  problem  just  discussed.  A two- 
phased  approach  was  employed.  In  the  first  phase,  the 
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question  of  what  probability  distribution  fits  the  failure 
data  was  to  be  settled.  If  the  Poisson  was  discredited,  a 
second  phase  was  to  be  undertaken  to  determine  whether 
incorporation  of  an  alternative  distribution  significantly 
changes  the  computed  requirement.  The  final  aspect  of 
the  problem  would  have  required  the  collection  and  evalu- 
ation of  so  much  additional  failure  data  that  it  appeared 
reasonable  not  to  do  so  unless  there  were  indications  from 
this  research  to  warrant  the  effort. 


Justification 


It  can  be  assumed  that  a requirements  computation 
reflecting  the  actual  demand  distribution  exprienced  by 
an  item  would  be  more  accurate  than  one  reflecting  an 
inappropriate  distribution.  A more  accurate  computation 
would  have  one  of  two  results.  In  the  first  case,  a 
smaller  requirement  would  be  computed,  saving  the  Air 
Force  the  money  it  would  have  spent  for  repair  or  pro- 
curement. In  the  second,  a larger  requirement  would  be 
computed,  necessitating  more  repair  and  procurement  funds. 

Improved  accuracy,  while  desirable,  is  not  free. 
In  order  to  implement  this  kind  of  improvement,  it  is 
possible  that  the  Air  Force  would  have  to  change  the  data 
collection  system  that  feeds  the  requirements  process. 

The  current  system  collects  only  the  total  number  of 
failures  and  total  operating  time  for  a component 
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population  (8:6).  These  data,  as  previously  noted,  are 
used  in  determining  the  demand  rate.  However,  they  are  not 
sufficient  for  determining  the  distribution  from  which 
the  demands  are  taken  unless  that  distribution  is  defined 
by  only  one  parameter.  For  distributions  that  are  defined 
by  two  or  more  parameters,  unaggregated  data  would  have 
to  be  collected  in  such  a way  that  the  parameters  couid 
be  estimated  directly  from  the  data. 

The  costs  associated  with  providing  the  additional 
data  to  the  requirements  system  are  not  the  only  ones  that 
might  be  incurred.  The  cost  of  performing  the  commutation 
itself  could  be  increased.  A 1970  AFLC  study  cited  the 
computational  efficiency  of  using  the  Poisson  rather  than 
other  discrete  probability  distributions  * noted  that 
corrections  would  be  required  to  use  continuous  distribu- 
tions to  approximate  an  inherently  discrete  process  such 
as  demand  generation  (22:29).  Additional  computational 
costs  could,  therefore,  be  anticipated  if  any  probability 
distribution  other  than  the  Poisson  were  incorporated  into 
the  requirements  computation  algorithm. 

If  the  requirements  computation  were  not  sensitive 
to  the  differences  in  demand  embodied  in  alternative  proba- 
bility distributions,  significantly  improved  accuracy  would 
not  be  expected.  If  the  requirements  computation  were, 
however,  sensitive  to  these  differences,  the  magnitude  of 
the  change  would  influence  whether  an  improved  technique 
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would  pay  for  itself  in  savings  or  whether  its  effect 
would  warrant  the  costs  associated  with  its  implementation. 
The  Air  Force  Logistics  Management  Center,' Gunter  Air 
Force  Station,  Alabama,  has  requested  that  this  sensi- 
tivity be  investigated  (23)  . 

Objectives  of  the  Research 

This  research  had  three  objectives.  Specifically, 
the  research  sought  to: 

1.  Identify  the  distribution  that  best  describes 
the  occurrence  of  failures  for  each  IMU  within  a time  frame 
appropriate  to  the  requirements  computation. 

2.  Determine  how  to  modify  the  D041  computation 
to  incorporate  the  identified  distribution. 

3.  Determine  whether  che  computation  is  sensitive 
to  the  proposed  modification  through  a comparison  of 
requirements  for  three  inertial  measurement  units  computed 
using  both  the  standard  and  the  modified  D041  process. 

Research  Hypotheses 

1.  For  each  of  the  inertial  measurement  units 
under  study,  the  daily  occurrences  of  failure  do  not  fol- 
low a Poisson  distribution. 

2.  Spares  computations  are  sensitive  to  the 
probability  distribution  that  is  assumed  to  describe  the 
underlying  equipment  failure  pattern. 
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CHAPTER  II 


METHODOLOGY 

This  chapter  presents  the  methodology  used  for 
identification  of  the  distribution  that  best  describes 
the  occurrence  of  failures  for  each  of  the  three  IMUs. 

The  first  subject  to  be  addressed  is  the  organization  of 
the  data  which  was  used.  The  data  collection  plan,  includ 
ing  a discussion  of  the  G078C  data  processing  system  from 
which  the  data  were  taken,  will  follow.  Finally,  there 
will  be  a general  discussion  of  the  analysis  techniques 
employed . 

lopulation  Defined 

Failure  data  for  all  aircraft  . nertial  nv.asure- 
mer.t  units  (IMUs)  in  the  Air  Force  inventory  constitute 
the  universe  for  this  research.  Within  the  universe, 
three  discrete  populations  were  defined  as  failure  data 
applicable  to  the  FLIP  unit  used  on  the  C-5A,  the  LN-15 
unit  used  on  the  B-52G/H,  and  the  KT-73  unit  used  on  the 
A-7D/E  and  the  AC-130H. 

Failure  data  associated  with  each  IMU  population 
were  arranged  in  chronological  order  by  failure  date. 
Population  subsets  were  defined  by  partitioning  the 
data  by  base  period.  A base  period  was  defined  as  eight 
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consecutive  calendar  quarters.  The  eight  quarter  block 
was  chosen  as  the  basis  for  investigation  because  the 
requirements  computation  system  (D041)  uses  eight  quarters 
of  failure  history  for  its  demand  forecasts  (2:p.l-4) . 
Figure  1 illustrates  the  partition  of  failure  data  by  base 
period.  Base  Period  1 includes  the  eight  quarters  begin- 
ning with  the  first  one  in  which  an  IMU  failure  occurred 
(e.g,  in  Figure  1,  quarters  70-1  through  71-4).  In  Base 
Period  2,  the  first  quarter  of  the  previous  base  period 
is  dropped  and  another  quarter  added  at  the  end  (e.g., 
quarters  70-2  through  72-1) . This  process  is  repeated 
until  the  last  complete  eight-quarter  base  period  has 
been  partitioned. 

Census  data  on  the  occurrences  of  failure  for 
each  IMU  from  its  introduction  into  the  inventory  until 
30  June  1977  were  obtained.  Reasons  for  establishing  the 
30  June  cutoff  date  will  be  discussed  within  the  context 
of  the  data  collection  plan.  Failures  occurring  before 
the  first  quarter  in  which  at  least  twenty  failures 
are  recorded  were  eliminated  from  further  consideration. 
This  criterion  was  arbitrarily  established  to  minimize 
the  effects  of  system  startup  on  later  analysis  and  to 
delimit  the  number  of  base  periods  to  be  investigated. 
Figures  2,  3,  and  4 identify  the  data  that  was  used  in  the 
analysis. 
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FLIP 

LN-15 

KT-73 

1 July  1971 
through 

30  June  1977 

1 January  1974 
through 

30  June  1977 

1 January  1972 
through 

30  June  1977 

Fig.  2.  Inclusive  Dates  of  Data  for  Each  IMU 


FLIP 

LN-15 

KT-73 

1,108 

723 

2,065 

Fig.  3.  Number  of  Failures  in  Data  Base  for 

Each  IMU 


FLIP 

LN-15 

KT-73 

16 

6 

14 

Number  of  Complete  Base  Periods  Available 
for  Analysis  for  Each  IMU 


Data  Collection  Plan 


Unaggregated  failure  data  for  the  three  IMUs  are 
available  in  the  automated  Inertial  Guidance  Reliability 
and  Configuration  Program  for  Aircraft,  designated  the 
G078C  data  system.  The  G078C  data  base  contains  five  types 
of  records:  (1)  inertial  system  history,  (2)  test  data, 

(3  and  4)  subindenture  records,  and  (5)  nonserialized 
piece  parts  (4:p.2-l).  This  research  is  concerned  with 
the  data  contained  on  the  first  type  of  record  which  is 
established  each  time  a failed  IMU  is  returned  to  the 
depot  for  repair.  The  data  in  G078C  were  verified  by 
personnel  of  the  Aerospace  Guidance  and  Metrology  Center 
(AGMC) , Newark  Air  Force  Station,  Ohio,  and  are  considered 
to  be  accurate  by  engineers  of  the  Reliability  and  Opera- 
tions Division  (AGMC/SNO)  who  routinely  use  the  data 
(14;  19)  . 

Among  the  data  recorded  when  a type  1 record  is 
established  is  the  date  of  the  unit's  arrival  at  the 
depot.  Until  approximately  1 July  1977,  all  three  IMUs 
were  managed  under  SCARS  (Serialized  Control  and  Reporting 
System)  procedures.  As  long  as  SCARS  procedures  were 
in  effect,  the  date  of  arrival  reflected  in  the  G078C 
data  base  represented  the  actual  date  on  which  the  item 
arrived  at  the  depot  from  the  base  at  which  the  failure 
actually  occurred.  Since  that  time,  the  date  of  arrival 
recorded  in  the  G078C  has  reflected  the  date  on  which  the 
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unit  was  input  to  the  depot  work  center  for  repair  (14) . 

The  change  in  reporting  procedures  provided  a logical 
point  at  which  to  establish  the  data  cutoff. 

It  should  be  noted  that  the  date  of  arrival  at 
the  repair  depot  is  not  the  same  as  the  date  of  failure. 

The  possibility  of  directly  obtaining  the  failure  date 
was  considered.  The  date  on  which  the  failed  IMU  was 
removed  from  an  aircraft  is  reported  on  AFTO  Form  88, 
Historical  Record  for  Inertial  Navigation  Unit  (24:p.2-4). 

A complete  file  of  AFTO  Forms  38  is  maintained  at  AGMC 
by  the  Reliability  and  Operations  Division  (14).  It 
was  determined  that  the  effort  required  to  extract  the 
data  from  the  forms,  validate  them  against  the  G078C 
data,  and  manually  build  new  data  files  would  he  excessive. 

As  an  alternative  to  obtaining  the  date  of  fail- 
ure diroctly  from  AFTO  Forms  88 . the  possibility  of  adjust- 
ing the  arrival  date,  as  recorded  in  GG78C,  to  approxi- 
mate the  actual  failure  date  was  also  considered.  The 
SCARS  Management  Indicator  Roport  (RCS;LOG-LOM(/)  71224) , 
prepared  quarterly  by  the  SCARS  monitor  at  the  Air  Logis- 
tics Center  at  which  the  IMU  is  managed,  reflects  the  aver- 
age number  of  days  spont  in  transit  from  base  to  depot 
for  each  IMU  (26;p.l3-6) . Adjusting  all  arrival  dates 
within  a quarter  by  the  same  amount  would  only  have 
shifted  the  data  and  would  not  have  yielded  useful 


information  concerning  actual  failure  dates.  This  alterna- 
tive, therefore,  was  rejected. 

Due  to  the  fact  that  the  actual  failure  dates  are 
not  readily  available,  it  was  assumed  that  the  pattern  of 
unserviceable  arrivals  at  the  depot  accurately  reflects 
the  distribution  of  failures  as  they  actually  occurred. 

Data  Preparation 

Arrival  dates  were  extracted  from  the  G078C 
type  1 records.  These  dates  are  located  in  columns  17 
through  21  or  the  standard,  eighty-column,  single-card 
record  (4:A5-1)„  Three  files  of  data  were  established, 
one  for  each  I MU . Each  of  these  files  was  sorted  chrono- 
logically in  order  to  partition  the  data  into  base  periods. 
Once  sorted  in  this  manner,  multiple  arrivals  on  a given 
date  were  indicated  by  multiple  occurrences  of  that  date 
on  the  file.  Individual  files  for  each  base  period  were 
extracted  from  the  three  chronological  master  files. 

Each  of  the  base  period  files  was  then  processed 
using  a computer  program  developed  by  the  authors.  The 
program  accomplished  the  following: 

1.  The  frequency  with  Which  each  date  appeared 
on  the  file  was  recorded. 

2.  Dates  which  did  not  appear  on  the  file  were 
recorded  with  a frequency  of  zero. 
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3.  All  dates  occurring  on  Saturday  and  Sunday 
were  eliminated  from  the  file.1 

The  resulting  data  files  were  in  the  format  required  for 
the  analysis  procedure  described  in  the  following  section. 

Research  Design 

The  hypothesis  that  the  number  of  arrivals  per 

day  in  each  base  period  follows  a Poisson  distribution 

was  tested  using  SIMFIT,  a computer  program  on  the 
2 

CREATE  system.  SIMFIT,  written  in  FORTRAN  IV,  uses  a 
histogram  construction  technique  to  compare  theoretical 

3 

distributions  to  data  supplied  by  the  user. 

Two  nonparametric  Goodness-of-Fit  tests,  the 
Kolmogorov-Smirnov  One  Sample  (K-S)  and  the  Chi-Square 
(X2 3),  are  employed  in  SIMFIT  to  determine  whether  a given 
theoretical  distribution  adequately  describes  the  user's 


1With  the  exception  of  four  arrivals  which  occur- 
red on  weekends  in  the  KT-73  master  file,  no  arrivals  were 
found  on  either  Saturday  or  Sunday  in  any  of  the  data  files. 
Since  arrivals  at  the  depot  could  not  bo  expected  on  week' 
ends,  it  did  not  appear  reasonable  to  include  weekend  dates 
since  they  accounted  for  over  25  percent  of  tho  dates  in 
a base  period. 

2 

CREATE  is  an  acronym  for  Computer  Resources  for 
Engineering,  Analysis,  Training,  and  Education.  The  sys- 
tem is  located  at  Wright-Patterson  AFD  OH  and  supports 
AFLC  and  Air  Force  Institute  of  Technology  (AFIT)  activi- 
ties . 


3 

Appendix  A provides  a brief  description  of  tho 
SIMFIT  computer  program. 
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data.  These  two  particular  tests  were  selected  for  inclu 
sion  in  the  program  because  of  their  distribution-free 
property.  The  K-S  test  is  preferred  over  the  x2  test 
because  of  its  relatively  greater  power  in  determining 
goodness-of-fit  (25:Ch.2) . One  limiting  factor  of  the 
K-S  test  is  that  it  is  not  applicable  when  the  population 
parameters  cannot  be  specified  in  advance  of  the  test. 

The  x2  test  should  normally  be  used  when  the  population 
parameters  must  be  estimated  from  sample  data  (12:86). 

SIMFIT  determines  the  mean  and  variance  of  the 
input  data  and  uses  them  to  calculate  the  parameter (s) 
of  the  distribution  being  tested  unless  the  user  speci- 
fies the  parameter (s) . In  this  research,  the  data  for 
each  base  period  represented  a subset  of  the  population 
rather  than  a sample  from  the  population.  Therefore,  the 
mean  and  variance  of  the  data  were  considered  to  be  popu- 
lation parameters  rather  than  estimates  derived  from  a 
sample,  and  it  was  considered  appropriate  to  use  the  K-S 
test. 

In  addition,  SIMFIT  assumes,  for  the  purposes  of 
the  x2  test,  the  parameters  are  always  estimated  from  a 
sample  and  reduces  the  degrees  of  freedom  for  the  test 
statistics  commensurate ly  (1) . For  these  reasons,  the 
K-S  test  was  used  exclusively  as  the  measure  of  the 
goodness-of-fit  for  this  research. 


The  results  of  the  K-S  test  were  used  to  accept 
or  reject  a particular  distribution  as  the  underlying 
distribution  for  the  arrivals  in  a given  base  period. 
SIMFIT  offers  a choice  of  three  levels  of  confidence 
for  the  G-O-F  tests:  90  percent,  95  percent,  and  99  per- 
cent (25:2).  A level  of  90  percent  was  chosen  as  appro- 
priate because  it  is  the  most  conservative  of  the  three. 
This  means  that  &=.10  offers  the  greatest  opportunity 
available  to  reject  a particular  distribution  for  a 
given  set  of  data . 

The  statistical  hypothesis  tested  is  stated  as 

follows: 


Null  Hq:  x ^ the  hypothesized  distribution 

with  the  appropriate  param- 
eter (s)  . 

Alternate  x $ the  hypothesized  distribution 

with  the  appropriate  param- 
eter (s)  . 

where  x is  the  number  of  arrivals  per  day.  The  null 
hypothesis  was  rejected  at  the  90  percent  level  of  confi- 
dence using  the  K-S  G-O-F  test  if  the  calculated  value  of 
D (the  K-S  statistic)  was  greater  than  the  critical  value 
of  0..  The  criterion  for  identifying  the  theoretical 
distribution  that  best  describes  the  distribution  of 
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arrivals  for  each  IMU  population  was  established  as  the 
distribution  that  passed  the  K-S  test  for  the  greatest 
number  of  base  periods. 

Research  Assumptions 

The  following  assumptions  apply  to  this  research: 

1.  Individual  IMU  failures  are  independent,  but 
the  rate  at  which  failures  occur  may  change  over  time. 

2.  The  data  extracted  from  the  G078C  data  base 
are  accurate. 

3.  SIMFIT  reliably  fits  known  distributions  to 
sample  data  and  accurately  estimates  the  parameters  of 
those  distributions. 

4 . The  distribution  of  unserviceable  IMU  arri- 
vals at  the  depot  accurately  reflects  the  occurrence  of 
IMU  failures. 

5.  The  subset  of  a population  may  be  treated 
as  a discrete  population  for  statistical  purposes. 

Research  Limitations 

The  following  limitations  to  the  research  effort 
have  been  identified: 

1.  There  has  been  no  attempt  to  establish  that 
the  three  IMUs  studied  are  representative  of  all  aircraft 
IMUs,  therefore,  the  findings  of  this  research  should  not 
be  generalized  without  further  corroboration. 
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2.  Findings  are  limited  to  an  indication  of 
whether  there  is  a problem,  in  terms  of  the  accuracy  of 
the  computation,  associated  with  the  using  of  the  Pois- 
son instead  of  a more  appropriate  distribution.  No 
attempt  has  been  made  to  assess  the  number  of  items 
involved,  nor  to  determine  the  costs  associated  with 
rectifying  any  problem.. 


CHAPTER  III 


ANALYSIS  AND  FINDINGS 


The  methodology  for  determining  the  distribution 
that  best  fits  the  failure  data  included  in  each  of  the 
base  period  data  files  was  discussed  in  the  last  chapter. 
The  approach  taken  was  first  to  determine  whether,  in 
fact,  the  Poisson  distribution  passed  the  K-S  test  for 
any  of  the  base  periods.  Once  that  was  done,  other  dis- 
tributions that  might  describe  the  data  were  examined 
using  the  SIMFIT  computer  program. 


The  Poisson  Distribution 


The  first  hypothesis  addressed  was  that  the 
occurrences  of  failure  do  not  follow  a Poisson  distribution 
This  hypothesis  was  first  examined  with  respect  to  the 
FLIP  IMU,  A complete  analysis  of  each  of  the  sixteen 
base  periods  revealed  that  none  of  these  subsets  of  data 
could  be  described  as  the  Poisson  distribution.  In  most 
cases,  the  difference  between  the  computed  value  of  the 
K-S  test  statistic  (MAX  D)  and  the  critical  value 
(Dcrit)  was  substantial.  To  illustrate  the  general 
lack  of  fit,  the  values  of  the  test  statistics  for  the 
tests  on  the  Poisson  distribution  for  the  FLIP  IMU  are 
summarized  in  Table  1.  The  fact  that  the  data  for  each 


TABLE  1 


VALUES  OF  THE  TEST  STATISTIC  COMPUTED  FOR  THE  SIMFIT 
K-S  TEST  IN  THE  POISSON  DISTRIBUTION  ANALYSIS 


<Dcrit  ‘ °-053 

for  a = 0.10) 

Base 

Computed  Value  of  Test 

Statistic 

Period 

FLIP 

LN-15 

KT-73 

71-3 

0.064 

• • * 

• • • 

71-4 

0.070 

• » • 

• • • 

72-1 

0.081 

* • • 

0.222 

72-2 

0.075 

• • • 

0.206 

72-3 

0.076 

• • • 

0.208 

72-4 

0.082 

• • • 

0.227 

73-1 

0,090 

• • • 

0.224 

73-2 

0.103 

• • • 

0.216 

73-3 

0.113 

• • • 

0.216 

73-4 

0.110 

• • • 

0.215 

74-1 

0.120 

0.106 

0.217 

74-2 

0.123 

0.114 

0.225 

74-3 

0.131 

0.121 

0.237 

74-4 

0.139 

0.133 

0.238 

75-1 

0.150 

0.151 

0.246 

0.1 
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5-2 


L56 


0.157 


0.244 


base  period  failed  the  G-O-F  test  and  that,  with  only  a 
few  exceptions,  the  failure  margin  was  substantial,  indi- 
cated that  the  Poisson  distribution  could  not  be  accepted 
as  a reasonable  descriptor  for  any  of  the  FLIP  failure 
data. 

The  test  statistics  resulting  from  the  SIMFIT 
analyses  of  the  LN-15  and  KT-73  data  are  also  shown  in 
Table  1.  As  can  be  readily  observed  in  all  cases,  the 
differences  between  the  computed  value  of  the  test  sta- 
tistic and  the  critical  value  are  relatively  large,  indi- 
cating a very  poor  fit  of  the  Poisson  distribution  to  the 
base  period  data.  The  cumulative  result  of  these  thirty- 
six  rejections  of  the  statistical  hypotheses  (at  ct=0.10) 
was  that  the  first  research  hypothesis  could  be  conclusively 
supported.  That  is,  for  these  three  inertial  measure- 
ment units,  it  can  be  concluded  that  the  daily  occur- 
rences of  failure  do  not  follow  a Poisson  distribution. 

The  Negative  Binomial 
Distribution 

At  this  point  in  the  research,  the  quest  to  find 
a more  likely  candidate  for  the  underlying  failure  dis- 
tribution was  begun.  The  previous  SIMFIT  products  were 
helpful  in  this  effort.  Among  the  statistics  computed 
in  the  program  for  each  set  of  data  were  the  arithmetic 
mean,  the  variance,  and  the  variance- to-mean  (V/M)  ratio. 


The  V/M  ratio  was  the  key  to  finding  a suitable  distribu- 
tion to  describe  the  IMU  failure  data. 

Three  theoretical  distributions  are  loosely  related 
through  the  V/M  ratio.  These  are  the  Poisson,  the  binomial, 
and  the  negative  binomial.  The  latter  two  are  members  of 
the  Bernoulli  family  of  distributions.  Considering  only 
these  three  distributions,  for  the  moment,  the  following 
properties  can  be  used  as  a guide  in  choosing  between 
them  (13:95)  i 

Negative  binomial  V/M  ratio  > 1 

Binomial  V/M  ratio  < 1 

Poisson  V/M  ratio  = 1 

Tables  2,  3,  and  4 show,  by  base  period,  the  mean, 

variance,  and  V/M  ratio  of  the  daily  failures  for  each  IMU 
as  computed  by  SIMFIT. 

The  V/M  ratio  for  each  set  of  base  period  data 
was  consistently  greater  than  one.  The  entire  data  base 
was,  therefore,  analyzed  by  SIMFIT  using  the  negative 
binomial  as  the  hypothesized  distribution.  The  results 
of  this  effort  are  shown  in  Tables  5,  6,  and  7 for  each 
IMU  by  base  period.  Of  the  thirty-six  data  sets  examined 
(sixteen  for  the  FLIP  IMU,  six  for  the  LN-15  IMU,  and 
fourteen  for  the  KT-73  IMU) , twenty-two  were  successfully 
fitted  with  a negative  binomial  distribution  at  a=0.10. 

For  each  base  period,  the  SIMFIT  program  computed 
the  arithmetic  mean  and  variance  data.  Using  these 
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TABLE  2 

DESCRIPTIVE  STATISTICS  BY  BASE  PERIOD 
FLIP  FAILURE  DATA 


Base 

Period 

Mean 

Variance 

V/M  Ratio 

71-3 

0.548 

0.755 

. 1.378 

71-4 

0.566 

0.830 

1.466 

72-1 

0.643 

1.000 

1.555 

72-2 

0.675 

1.014 

1.502 

72-3 

0.688 

1.051 

1.527 

72-4 

0.710 

1.107 

1.559 

73-1 

0.757 

1.190 

1.572 

73-2 

0.791 

1.357 

1.716 

73-3 

0.829 

1.510 

1.822  , 

73-4 

0.849 

1.603 

1.888 

74-1 

0.839 

1.595 

1.901 

74-2 

0.839 

1.621 

1.932 

74-3 

0.837 

1.656 

1.979 

74-4 

0 .834 

1.719 

2.061 

75-1 

0.807 

1.761 

2.182 

j 

i 


75-2 


0.793 


1.788 


2.254 


TABLE  3 


DESCRIPTIVE  STATISTICS  BY  BASE  PERIOD 
LN-15  FAILURE  DATA 


Base 

Period 

Mean 

Variance 

V/M  Ratio 

74-1 

0.692 

1.212 

1.752 

74-2 

0.757 

1.327 

,1.753 

74-3 

0.797 

1.404 

1.762 

74-4 

0.864 

1.578 

1.826 

75-1 

0.956 

1.924 

2.012 

75-2 

1.029 

2.076 

2.018 

TABLE  4 

DESCRIPTIVE  STATISTICS  BY  BASE  PERIOD 

KT-73  FAILURE  DATA 

Base 

Period 

Mean  Variance  V/M  Ratio 

72-1 

1.378 

4.052 

2.941 

72-2 

1.477 

4.174 

2.826 

72-3 

1.671 

4.670 

2.795 

72-4 

1.753 

5.036 

2.864 

73-1 

1.743 

4,748 

2.724 

73-2 

1.735 

4.580 

2.640 

73-3 

1.647 

4.277 

2.5P7 

73-4 

1.636 

C.227 

2.584 

74-1 

1.586 

3.340 

2.484 

74-2 

1.595 

4 . 01 2 

2.515 

74-3 

1.558 

3.912 

2.511 

74-4 

1.501 

3.744 

2.494 

75-1 

1.501 

3.889 

2.591 

75-2 

1.478 

3.748 

2.536 

i 

'i 
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TABLE  5 


NEGATIVE  BINOMIAL  SIMFXT  RESULTS  FOR  THE 
FLIP  FAILURE  DATA  BY  BASE  PERIOD 

(Dcrit  ~ 0.053  for  a *>  0.10) 


Base 

Period 

Max  D 

K-S 

Tost 

Negative  Binomial 
Parameters 

m 

P 

71-3 

0.004 

Pailed 

1 

0.726 

71-4 

0.044 

Passod 

1 

0.602 

72-1 

0.037 

Passed 

1 

0.643 

72-2 

0.001 

Failed 

1 

0.666 

72-3 

0.076 

Failed 

1 

0.655 

72-4 

0.060 

Failed 

1 

0.642 

73-1 

0.076 

Failod 

1 

0.636 

73-2 

0.026 

Passed 

1 

0.502 

73-3 

0.010 

Passod 

1 

0.549 

73-4 

0.026 

Passed 

1 

0.529 

74-1 

0.025 

Passod 

1 

0.526 

74-2 

0.037 

Passed 

1 

0.518 

74-3 

0.059 

Failed 

1 

0.505 

74-4 

0.008 

Failod 

1 

0.405 

75-1 

0.138 

Failod 

1 

0.458 

75-2 

0.164 

Failod 

1 

0.444 
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TABLE  6 


NEGATIVE  BINOMIAL  SIMFIT  RESULTS  FOR  THE 
LN-15  FAILURE  DATA  BY  BASE  PERIOD 

(D  ..  = 0.053  for  a = 0.10) 


3ase 

Period 

Max  D 

K-S 

Test 

Negative  Binomial 
Parameters 

m 

P 

74-1 

0.037 

Passed 

1 

0.571 

74-2 

0.017 

Passed 

1 

0.571 

74-3 

0.028 

Passed 

1 

0.568 

74-4 

0.027 

Passed 

1 

0.548 

75-1 

0.038 

Passed 

1 

0.497 

75-2 

0.022 

Passed 

1 

0.495 
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TABLE  7 


NEGATIVE  BINOMIAL  SIMFIT  RESULTS  FOR  THE 
KT-73  FAILURE  DATA  BY  BASE  PERIOD 


.34 

Failed 

81 

Failed 

38 

Passed 

50 

Passed 

31 

Passed 

30 

Passed 

28 

Passed 

29 

Passed 

43 

Passed 

43 

Passed 

49 

Passed 

60 

Failed 

83 

Failed 

78 

Failed 

s 
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statistics,  the  parameters,  m and  p,  of  a negative  binomial 
distribution  were  computed.  The  theoretical  distribution, 
thus  described,  was  then  compared  to  the  data  and  analyzed 
using  the  K-S  G-O-F  test.  The  parameters  of  the  distri- 
bution used  in  each  test  are  also  shown  in  the  tables  cited 
above . 


Complete  negative  binomial  test  results,  a detailed 
numerical  description  of  the  data,  and  the  associated  data 
histogram  for  each  base  period  for  each  IMU,  are  located 
in  appendixes  as  follow i 

Appendix  IMU 

B FLIP 

C LN-15 

D KT-73 


The  tables  and  histogram  presented  in  these 
appendixes  are  products  of  the  SIMFIT  computer  program. 
Note  that  some  of  the  tables  in  Appendix  D,  the  KT-73 
analyses,  have  the  tail  the  distribution  truncated  such 
that  the  cumulative  probability  does  not  reach  a value  of 
1.000.  This  was  done  to  facilitate  the  placement  of  the 
tables  on  the  page  and  does  not,  in  any  way,  affect  the 
test  results. 
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Other  Distributions 

Having  found  that  the  negative  binomial  distribu- 
tion adequately  deserved  a majority  (61.1  percent)  of 
the  base  period  data  sets,  the  research  effort  was 
extenucd  to  the  investigation  of  other  theoretical  distri- 
butions. Since  many  of  the  common  statistical  distribu- 
tions are  contained  in  the  SIMFIT  computer  program,  it  was 
relatively  easy  to  compare  each  of  them  to  any  given  set  of 
data. 

Rather  than  analyzing  every  data  set,  which  would 
have  required  a large  amount  of  computer  time,  it  was 
decided  to  sample  from  the  base  periods  for  each  IMU 
population.  Any  distributions  which  appeared  to  closely 
describe  the  data  could  then  be  investigated  further. 

The  base  periods  for  examination  were  randomly  selected 
for  each  IMU.  Since  the  FLIP  and  KT-73  IMUs  had  a rela- 
tively large  number  of  base  periods  from  which  selection 
could  be  made,  three  were  chosen  from  each  of  these 
populations.  Only  two  base  periods  were  chosen  for  the 
LN-15  IMU. 

Tables  8,  9,  and  10  contain  the  K-S  test 
results  (MAX  D values)  from  these  analyses.  Note  that 
none  of  the  test  statistics  except  some  of  those  for  the 
negative  binomial  runs,  are  less  than  the  critical  D value 
(Dcrit  = 0.053).  The  negative  binomial  K-S  test  results 
which  correspond  exactly  tc  those  presented  earlier,  are 
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TABLE  8 


K-S  TEST  STATISTICS  FOR  ALL  DISTRIBUTIONS  AVAILABLE 
IN  SIMFIT  FOR  THE  FLIP  FAILURE  DATA 

{Dcrit  = 0.053  for  a = 0.10) 


Distribution 

Test  Statistics 

FL-72-2 

FL-73-2 

FL-74-3 

Erlang 

• • • • 

• i • • 

• « • * 

Weibull 

0.412 

0.439 

0.433 

Gamma 

0.414 

0.441 

0.434 

Pearson  XI 

0.220 

0.204 

0.183 

Lognormal 

0.296 

0.263 

0.268 

Normal 

• • • • 

0.127 

0.170 

Uniform 

0.306 

0.313 

0.345 

Beta 

• • • • 

• • • • 

• • • • 

Triangular 

0.601 

0.573 

0.561 

Poisson 

0.075 

0.103 

0.131 

Negative  Binomial 

0.081 

0.026 

0.059 

Positive  Binomial 

• • • • 

• * * • 

• • • • 

32 


TABLE  9 


K-S  TEST  STATISTICS  FOR  ALL  DISTRIBUTIONS  AVAILABLE 
IN  SIMFIT  FOR  THE  LN-15  FAILURE  DATA 

^Dcrit  = 0.053  for  a = 0.10) 


Distribution 

Test 

Statistics 

LN-74-2 

LN-74-4 

Erlang 

• • • • 

• • • • 

Weibull 

0.414 

0.440 

Gamma 

0.415 

0.443 

Pearson  XI 

0.192 

0.178 

Lognormal 

0.298 

0.273 

Normal 

0.142 

0.160 

Uniform 

0.332 

0.325 

Beta 

• • • • 

• • • • 

Triangular 

0.573 

0.546 

Poisson 

0.114 

0.133 

Negative  Binomial 

0.017 

0.027 

Positive  Binomial 

• • • • 

• • • • 
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TABLE  10 


K-S  TEST  STATISTICS  FOR  ALL  DISTRIBUTIONS  AVAILABLE 


IN  SIMFIT 

FOR  THE  KT-73 

FAILURE  DATA 

(Dcrit 

= 0.053  for  a 

= 0.10) 

Distribution 

Test  Statistics 

KT-72-2 

KT-73-2 

KT-75-2 

Erlang 

0.088 

0.097 

0.113 

Wsibull 

0.548 

0.543 

0.498 

Gamma 

0.556 

0.552 

0.507 

Pearson  XI 

0.121 

0.127 

0.124 

Lognormal 

0.168 

0.152 

0.206 

Normal 

0.209 

0.289 

0.241 

Uniform 

0.376 

0.332 

0.348 

Beta 

• • • • 

• i • • 

• • • • 

Triangular 

0.459 

0.396 

0.414 

Poisson 

0.206 

0.211 

0.244 

Negative  Binomial 

0.081 

0.030 

0.078 

Positive  Binomial 

• • • * 

• • • • 

• • • • 
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shown  for  comparison  with  the  other  twelve  distributions. 
When,  in  the  tables,  no  value  is  indicated  for  a particu- 
lar distribution,  the  data  were  out  of  limits  for  the 
theoretical  distribution  and  no  SIMFIT  analysis  could 
be  performed.  On  the  basis  of  this  sample,  it  was  con- 
cluded that  none  of  the  distributions  examined,  except  the 
negative  binomial,  could  be  expected  to  describe  the  fail- 
ure data  for  these  IMUs.  Research  in  this  area  was,  then, 
discontinued. 

Summary  of  Findings 

This  phase  of  the  research  addressed  the  hypothesis 
that,  for  each  of  the  inertial  measurement  units  under 
study,  the  daily  occurrences  of  failure  do  not  follow  a 
Poisson  distribution.  In  attempting  to  fit  the  Poisson 
to  data  in  sixteen  base  periods  for  the  FLIP  unit,  six 
base  periods  for  the  LN-15  unit,  and  fourteen  base  periods 
for  the  KT-73  unit,  not  a single  case  was  identified  in 
which  that  distribution  could  be  accepted  as  a description 
of  the  data.  The  first  research  hypothesis  is  supported 
by  these  findings. 

The  negative  binomial  distribution  was  found  to 
fit  data  in  seven  of  sixteen  base  periods  for  the  FLIP 
unit,  all  six  base  periods  for  the  LN-15  unit,  and  nine 
of  fourteen  base  periods  for  the  KT-73  unit.  No  other 
distributions  were  found  to  fit  data  from  a sample  of 
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base  periods  for  each  IMU.  Therefore,  according  to  the 
criterion  established  for  this  research,  it  was  concluded 
that  the  negative  binomial  is  the  best  distribution  to 
describe  the  daily  occurrences  of  failure  for  each  of  the 
three  IMUs. 
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CHAPTER  IV 


APPLICATION  OP  FINDINGS 

The  second  phase  of  this  research  considered  the 
use  of  a negative  binomial  demand  distribution  in  the 
existing  requirements  computation  process.  The  last  two 
research  objectives  identified  in  Chapter  I were  addressed 
at  this  point.  The  first  objective  was  to  determine  how 
to  modify  the  D041  computation  to  incorporate  the  identi- 
fied demand  distribution.  The  other  was  to  determine 
whether  the  computation  would  be  sensitive  to  the  modifi- 
cation. 

This  chapter  begins  with  a review  of  some  of  the 
theory  behind  the  Variable  Safety  Level  (VSL)  feature  of 
the  D041  computation.  A discussion  of  the  complexities 
surrounding  incorporation  of  a negative  binomial  demand 
distribution  in  VSL  follows.  These  complexities  pre- 
cluded performance  of  a sensitivity  analysis  to  meet  the 
final  objective  of  the  research. 

Computation  of  Requirements 

The  Variable  Safety  Level  (VSL)  feature  of  the 
D041  requirements  computation  is  based  upon  the  METRIC 


model. 


METRIC  is  a model  for  determining  both  require- 
ments and  distribution  of  recoverable  items  in  a two- 
echelon  inventory  system.  The  objective  of  the  model 
is  to  determine  the  base  and  depot  stock  levels  which 
minimize  total  expected  backorders  for  a specific  set 
of  items  and  bases  subject  to  an  investment  constraint 
....  A backorder  occurs  at  a point  in  time  at  which 
there  is  an  unsatisfied  demand  at  base  level  [17:6]. 


Backorder  Computation.  Most  aircraft  recoverable 
spare  parts  are  included  in  a class  of  items  having  a 
high  unit  cost  or  a low  demand.  The  optimal  inventory 
policy  associated  with  this  class  of  items  is  "to  place 

a reorder  immediately  whenever  a demand  occurs  [11:1]." 

4 

With  this  policy,  the  state  of  the  inventory  at  a given 
time  is  defined  by  the  number  of  units  in  resupply.  Units 
that  have  been  reordered  (requisitioned  from  the  depot  or 
duo  in  from  base  maintenance)  are  considered  to  be  in 
resupply  until  delivered  to  base  supply  in  serviceable 
condition.  When  the  number  of  items  in  resupply  exceeds 
the  spare  (safety)  stock  that  is  in  the  system,  there 
are  backorders  (11:2). 

In  order  to  determine  the  level  of  spare  stock 
that  yields  the  minimum  expected  backorders,  the  steady 

5 

state  probabilities  for  the  number  of  units  in  resupply 

4 

State  is  variously  defined  as  a specific  measure- 
able  condition  of  a system  (c:314)  or  as  a specific  value 
for  the  random  variables  of  a stochastic  process  (6:582) . 

5 

The  concept  of  steady  state  comes  from  queueing 
theory,  "In  the  steady  state,  probabilities  are  inde- 
pendent of  time  [6:468]."  On  the  other  hand,  in  the  trans- 
ient state,  probabilities  are  dependent  upon  the  initial 
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must  be  used.  These  probabilities  are  also  known  as  the 
state  probabilities,  and  can  be  alternatively  expressed 
as  the  probabilities  of  the  number  of  demands  in  a time 
interval  of  specified  length.  Expected  backorders  are 
computed  by  the  formula: 


B (s) 


(x-s)h(x) 

\=s+l 


where : 

s is  the  stock  level 

h(a)  is  the  steady  state  probabilities  of  a:  units  in 
resupply. 

In  VSL,  h(x)  is  a negative  binomial  probability  distri- 

i 

but ion  (7).  ! 


Logarithmic  Poisson  Demand . The  literature 
identifies  two  different  demand  processes  that  have  nega- 
tive binomial  state  probabilities . The  more  thoroughly 
documented  of  the  two  processes  is  that  in  which  demand 
is  assumed  to  occur  according  to  a logarithmic  Poisson 
process.  The  logarithmic  is  one  of  a family  of  compound 

conditions.  For  an  introduction  to  queueing  theory, 
see  Principle s of  Operations  Research  for  Management , 
Chapter  12  (6) . 
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Poisson  processes  that  are  generalizations  of  the  simple 
Poisson  process  introduced  in  Chapter  I.  With  a compound 
Poisson  process,  customers  arrive  (i.e.,  place  demands), 
according  to  a (simple)  Poisson  process.  If  each  customer 
can  demand  more  than  one  unit  at  a rime,  the  number  of 
demands  occurring  within  a specified  interval  of  time  has 
a compound  Poisson  distribution.  If  customers  are  limited 
to  demanding  one  unit  at  a time,  demands  occur  according  to 
a simple  Poisson  distribution  (20:4).  In  either  case,  the 
distribution  of  time  between  events,  batches  of  demand  for 
the  compound  process  or  individual  demands  for  the  simple 
process,  is  exponential  (11:4). 

The  logarithmic  Poisson  distribution  is  obtained 
by  considering  batches  of  demand  where  the  number  of 
hatches  (customer  arrivals]  follows  a Poisson  process 
and  the  number  of  demands  has  a logarithmic  distri- 
bution [21:8]. 

The  RAND  literature  cites  a theorem  regarding 
system  performance  under  the  inventory  policy  for  recov- 
erable items: 

...  if  demand  is  Poisson,  then  the  number  of  , 
units  in  resupply  in  the  steady  state,  x,  is  also 
Poisson  for  any  distribution  of  resupply  [times] . 

The  Poisson  state  probability  depends  on  the  mean  of 
the  resupply  distribution,  but  not  on  the  distri- 
bution itself  [11:2]. 

For  a simple  Poisson  process,  with  mean  demand  rate  (mean 
number  of  demands  per  iterval  of  time),  X,  and  mean  resup- 
ply time,  T,  the  state  probabilities  are  simple  Poisson 
with  rate  XT  (11:3).  Similarly,  for  a compound  Poisson 
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process  with  ouotomsv  aw&oal  vato*  X,  and  moan  resupply 
time,  T,  the  stato  probabilities  aro  oompound  Poisson  with 
rate  \T  (11:7)  . Sherbrooke  demonstrated  that  the  atato 
probabilities  derived  from  a logarithmic  Poiaoon  domand 
proco33  can  bo  oxprossod  ao  a negative  binomial  distri- 
bution  (20:12). 

Dayosian  Inference.  In  the  same  discussion, 
Sherbrooke  acknowledged  that  "there  aro  other  non-compound 
Poiason  processes  that  aloe  yield  negative  binomial  atato 
probabilities  (20:12)."  One  such  process  is  "the  proba- 
bility distribution  of  domand,  whore  time  is  broken  into 
intervals  of  fixed  longth  and  domand  is  Poisson-distributed 
with  poriod  moans  selected  from  a gamma  distribution 
[20:12-13]."  This  statement  would  b^  of  only  passing 
interest  wore  it  not*  for  the  concept  of  Dayosian  inference 
that  is  prominent  in  the  litoraturo.  Four  RAND  documents, 
RM-1413  (27),  RM-43G2-PR  (9),  RM-4720-PR  (10),  and  RM-507B- 
PR  (21) , argue  that  Dayosian  inference  con  bo  applied  in 
demand  analysis . The  discussion  that  follows  describes 
this  application  of  Dayosian  inference. 

It  in  argued  that  observed  mean  domand  for  an 
itom,  x,  is  insufficient  to  determine  the  true  moan  demand 
for  t>H3  item,  0.  Using  the  Dayosian  approach,  "wo  can 
increase  our  knowledge  of  an  item  by  analysing  the  behavior 
of  related  items  in  tho  supply  system  (10:6)."  If  the 
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specified  item  is  considered  to  be  one  of  a large,  finite 
group  of  items  (i.e.,  the  recoverable  item  inventory), 
each  with  a true  mean  demand,  0,  there  is  a distribution  of 
0 values  for  the  whole  group.  The  distribution  of  0 
values  becomes  the  Bayesian  prior  probability  distribution 
(9:3).  Since  th$  true  mean  demands  for  the  individual 
items  are  not  known,  the  distribution  of  0 values  must 
be  estimated. 

Estimation  of  the  prior  distribution  requires  that 
the  form  of  the  distribution  be  specified  and  that  the 
distribution's  parameters  be  estimated.  The  form  of  the 
prior  distribution  is  chosen  to  reflect  what  is  known 
about  the  inventory  as  a whole — that  "most  items  have  low 
demand  while  a few  may  have  very  high  demand  [9:5]." 
Therefore  the  distribution  chosen  should  accommodate  any 
nonnegative  real  value  for  0 and  be  skewed  to  the  right. 

The  data  used  to  estimate  the  parameters  of  the  prior 
distribution  were  identified  as  a cross  section  of 
observed  demand.  This  set  of  data  includes  the  total 
demand  observed  over  a specified  time  period  for  each 
item  in  the  inventory.  It  is  desirable  that  the  esti- 
mated distribution  have  at  least  two  parameters  "so  that 
the  first  two  moments  of  this  prior  distribution  can  be 
estimated  from  the  cross-sectional  data  [9:5]." 
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The  essence  of  the  Bayesian  approach  is  that: 

. . . the  system  of  inventory  items  provides  the 
information  for  a prior  distribution  common  to  all 
items  [and  this  information  is  combined]  with  the 
demand  data  on  a specific  item  to  give  a posterior 
distribution  for  that  item  [10:7-8], 

Several  forms  for  the  prior  distribution  were  pro- 
posed in  the  literature,  however,  the  gamma  distribution 
was  specified  for  the  METRIC  model.  It  should  be  recog- 
nized that  when  the  prior  distribution  is  assumed  to  be 
gamma  and  the  demand  distribution  is  assumed  to  be  loga- 
rithmic Poisson,  the  posterior  probabilities  that  can  be 
computed  are  only  approximations.  However,  when  the  demand 
distribution  is  assumed  to  be  simple  Poisson,  the  posterior 
probabilities  can  be  computed  exactly  and  yield  a negative 
binomial  distribution  (21:31). 

I 

Implications  of  Research 
Findings 

The  fact  that  the  VSL  computation  incorporates 
a negative  binomial  probability  distribution  to  compute 
expected  backorders  could  be  a consequence  of  one  of 
two  lines  of  reasoning.  One  accepts  the  technique  of 
Bayesian  inference  and  assumes  the  demand  distribution 
to  be  simple  Poisson  and  the  prior  distribution  to  be 
gamma.  The  other  rejects  the  technique  of  Bayesian  infer- 
ence and  assumes  the  demand  distribution  to  be  logarithmic 


Poisson. 


The  research  methodology  specifically  addressed 
the  assumption  that  demand  for  three  IMUs  follows  a simple 
Poisson  distribution.  This  assumption  was  not  supported 
by  the  research.  It  has  been  suggested,  however,  that  the 
fact  that  the  data  were  aggregated  at  daily  intervals  may 
have  masked  the  existence  of  a simple  Poisson  demand 
process  (18) . 

The  assumption  of  a logarithmic  Poisson- demand 
process  was  not  addressed  because  SIMFIT  does  not  include 
a logarithmic  (or  any  other  compound)  Poisson  distribu- 
tion. Although  provisions  are  available  for  specifying 
additional  distributions,  doing  so  is  a major  program- 
ming and  validation  effort  which  could  not  be  completed 
within  the  time  constraints  imposed  on  the  research. 
However,  a comment  about  this  assumption  can  be  made. 
Sherbrooke  established  that  "any  compound  Poisson  dis- 
tribution with  a positive,  discrete  compounding  dis- 
tribution has  a variance  that  equals  or  exceeds  its  mean 
[20:5]."  It  will  be  recalled  that  the  demand  data  for  all 
of  the  base  periods  have  variance-to-mean  ratios  greater 
than  one  (i.e,,  the  variance  exceeds  the  mean).  While 
this  fact  does  not  necessarily  imply  that  the  demand 
distributions  are  compound  Poisson,  that  possibility  can- 
not be  dismissed. 

The  conclusion  that  the  demand  distribution  for 
the  three  IMUs  is  negative  binomial  presents  a problem. 
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The  literature  does  not  indicate  the  distribution  of  the 
state  probabilities  associated  with  negative  binomial 
demand  distributions.  It  was  beyond  the  present  capabili- 
ties of  the  authors  to  identify  the  appropriate  state 
probability  distribution.  Therefore,  no  analysis  could 
be  made  of  the  sensitivity  of  the  expected  backorder 
computation  to  a substitution  of  distributions. 
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CHAPTER  V 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  conclusions  drawn  from  the  research  on  the 
failure  data  for  the  FLIP,  LN-15,  and  KT-73  inertial  mea- 
surement units  are  discussed  in  this  final  chapter.  Each 
research  hypothesis  is  addressed  in  turn.  The  chapter 
concludes  with  several  recommendations  for  future  research. 


The  First  Research  Hypothesis 


The  hypothesis  first  addressed  m this  research 
was  that:  For  each  of  the  inertial  measurement  units  under 
study , the  daily  occurrences  of  failure  do  not  follow  a 
Poisson  distribution.  Statistical  tests  for  goodness-of- 
fit  of  the  Poisson  were  performed  on  failure  data  parti- 
tioned by  base  period.  The  Poisson  distribution  could  not 
be  fit  to  any  of  the  base  period  data  sets.  These  results 
strongly  support  this  hypothesis. 

Additional  investigation  indicated  that  the  nega- 
tive binomial  distribution  adequately  described  the  data 
in  61.1  percent  of  the  base  periods.  None  of  the  other 
distributions  available  in  SIMFIT  fit  any  of  the  data.  In 
consonance  with  the  criterion  stated  at  the  outset  of 
this  research  effort,  it  is  concluded  that  the  negative 


binomial  is  the  best  distribution  for  describing  the 
occurrences  of  failure  for  each  of  the  IMUs  studied. 

The  Second  Research  Hypothesis 

The  second  hypothesis  addressed  was:  Spares 
computations  are  sensitive  to  the  probability  distribution 
that  is  assumed  to  desox'ibe  the  under \ failure  pattern. 

No  conclusions  can  be  drawn  from  this  j:ec  ■jarch  that  either 
support  or  refute  this  hypothesis  because  i sensitivity 
analysis  could  not  be  performed. 

It  was  reported  that  the  objective  of  the  VSL 
computation  is  to  determine  the  base  and  depot  stock  levels 
that  minimize  total  expected  backorders,  subject  to  a 
monetary  constraint.  The  sensitivity  of  this  computation 
to  a change  in  the  demand  distribution  would  be  determined 
by  comparing  the  optimal  stock  levels  and  the  associated 
number  of  expected  backorders  that  result  from  using 
different  demand  distributions.  In  order  to  perform  the 
computations,  it  is  necessary  to  identify  the  steady 
state  probabilities  for  the  number  of  items  in  resupply 
for  the  demand  distributions  under  consideration.  A 
review  of  the  literature  failed  to  reveal  the  state  proba- 
bilities associated  with  a negative  binomial  demand  distri- 
bution. Without  this  critical  piece  of  information,  the 
analysis  was  impossible. 
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Recommendations 

It  was  concluded  that,  in  the  case  of  three  IMUs, 
it  is  inappropriate  to  assume  the  existence  of  a simple 
Poisson  demand  distribution.  Ideally,  this  research  would 
have  identified  a modification  to  the  D041  computation, 
specifically  to  the  VSL  segment,  that  incorporates  a more 
appropriate  demand  distribution,  the  negative  binomial. 
Without  that  modification,  and  subsequent  comparison  of 
its  results  with  those  based  on  the  present  assumption,  no 
recommendations  regarding  the  D041  can  bo  made.  However, 
the  research  findings  do  suggest  some  avenues  for  further 
research. 

1.  The  assumption  of  a logarithmic  Poisson  demand 
distribution  should  be  examined.  Srnce  the  failure  data 
in  all  base  periods  exhibited  variances  greater  than  means, 
one  characteristic  of  this  distribution  is  established. 

This  research  would  require  development  of  a logarithmic 
Poisson  subroutine  for  SIMFIT.  The  type  1 records  in  the 
G078C  data  system  identified  the  base  at  which  each  failure 
occurred,  so  identifying  the  distribution  of  customer 
arrivals  and  the  distribution  of  demands  per  customer  should 
not  be  difficult. 

2.  The  avenue  of  research  that  appears  to  be  most 
challenging  is  the  identification  of  the  state  probability 
distribution  associated  with  a negative  binomial  demand 
distribution.  It  will  be  recalled  that  queueing  theory 


provided  some  of  the  concepts  on  which  METRIC,  and  later 
VSL,  was  based.  Budnick  et  al.,  suggest  that  derivation 
of  analytic  solutions  to  queueing  problems  poses  enormous 
problems  to  all  but  the  most  accomplished  of  mathema- 
ticians. They  recommend  obtaining  the  solution  by  simu- 
lation (6:439). 

3.  Research  should  be  undertaken  to  determine 
whether  the  finding  that  demand  is  best  described  by  a 
negative  binomial  distribution  applies,  as  well,  to  other 
IMUs.  That  these  IMUs  are  installed  on  aircraft  with 
such  diverse  characteristics  suggests  that  the  finding  may 
be  generic  to  IMUs  rather  than  attributable  to  some 
peculiarity  of  their  operational  environments. 

4.  The  G078C  system  records  much  more  data  than 
was  used  in  this  research.  One  item  that  can  be  obtained 
from  the  data  base  is  the  number  of  hours  the  IMU  operated 
before  it  failed.  This  type  of  data  could  be  used  in 
making  an  actuarial  forecast  of  failures.  Since  the 
assumptions  of  the  D041  computation  seem  to  be  violated, 

it  may  be  more  reasonable  to  adopt  the  actuarial  technique, 
which  is  currently  used  by  the  Air  Force  only  to  develop 
aircraft  engine  requirements,  rather  than  attempting  to 
modify  the  D041.  The  efficacy  of  this  approach,  of  course, 
would  depend  on  whether  these  findings  apply  across-the- 
board  to  all  IMUs  or  only  to  these  three. 
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APPENDIX  A 
SIMFIT 
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SIMFIT  is  a computer  program,  written  in  FORTRAN  IV 
which  is  used  to  compare  data  to  theoretical  probability 
distributions.  Each  distribution  available  in  the  program 
is  handled  by  a separate  subroutine  which  operates  inde- 
pendently from  the  main  program  and  the  other  distri- 
bution subroutines.  The  following  distributions  are  avail 
able  in  the  SIMFIT  program  (1): 

Cumulative- 

Erlang 

Weibull 

Gamma 

Pearson  XI 

Lognormal 

Normal 

Uniform 

Beta 

Triangular 

Poisson 

Negative  Binomial 
Positive  Binomial 

In  addition  to  the  thirteen  distributions  listed  above, 
the  program  provides  for  the  inclusion  of  up  to  seven 
subroutines  which  must  be  written  by  the  user  to  specify 
other  theoretical  distributions  or  empirical  distributions 
of  his  own  design. 


Goodness-of-Fit  Tests 

The  program  employs  both  the  Kolmogorov-Smirnov 
(K-S)  and  the  Chi-Square  (x2)  tests  for  goodness-of-fit 
(G-O-F) . The  tests  are  used  to  determine  whether  the 
values  of  a sample  can  reasonably  be  thought  to  have 
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come  from  a population  described  by  the  selected  theo- 
retical distribution.  The  SIMFIT  user  should  familiarize 
himself  completely  with  the  mechanics  and  theory  behind 
each  of  these  G-O-F  tests  prior  to  performing  any  analysis 
with  SIMFIT. 

It  should  be  noted  that  SIMFIT  uses  each  test 
exactly  as  it  is  used  in  classical  nonparametric  statistic 
The  K-S  test  treats  each  observation  separately  and  does 
not  cause  a loss  of  information  by  combining  cells  which 
have  small  numbers  of  observations.  The  x2  test,  on  the 
other  hand,  combines  cells  which  contain  less  than  five 
observations  and  reduces  the  degrees  of  freedom  used  in 
determining  the  x2  statistic  accordingly. 

SIMFIT  has  statistical  tables  for  both  G-O-F 
tests  written  directly  into  the  program,  but  only  the 
values  for  three  confidence  levels  are  included:  90,  95, 
and  99  percent.  If  the  user  desires  to  use  a different 
confidence  level,  he  must  change  these  tables  accordingly. 
This  io  done  by  revising  the  appropriate  lines  of  coding 
to  reflect  the  desired  values  (1). 

Attractive  Features  of  SIMFIT 

Since  the  program  is  written  in  FORTRAN  IV,  it  can 
be  easily  extended  and  modified  to  fit  the  specific  needs 
of  the  user.  Care  must  be  taken,  however,  to  investigate 
fully  the  impact  of  a change  in  SIMFIT  on  other  parts  of 
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the  program.  For  this  reason,  the  user  should  become 
thoroughly  familiar  with  the  program  listing,  which  is 
reasonably  well  documented,  before  attempting  to  modify 
the  coding. 

The  output  is  formatted  in  such  a manner  that 
it  will  fit  easily  into  an  8-1/2  by  11  inch  format.  This 
feature  is  especially  useful  for  the  thorough  documen- 
tation of  a research  effort. 

SIMFIT  provides  an  intermediate  printout  of  some 
of  the  input  variables  as  well  as  some  statistical  informa- 
tion drawn  from  the  data.  Items  included  in  the  statistics 
are  the  mean,  standard  deviation,  variance,  variance-to- 
mean  ratio,  and  the  third  and  fourth  moments  about  the 
mean.  This  information  is  valuable  in  accurately  describ- 
ing and  identifying  the  data. 

As  a result  of  the  computer  analysis,  specific 
numerical  information  is  printed  out  for  each  cell  of  the 
data,  histogram.  The  user  will  find  this  information 
essential  for  detailed  data  analysis. 

In  addition  to  the  numerical  information  pro- 
vided, the  program  produces  a histogram  of  the  data  values. 
For  convenience  in  picturing  the  shape  of  the  theoretical 
distibution  being  investigated,  an  approximation  of  the 
distribution  is  superimposed  onto  the  data  histogram. 

While  the  whole  picture  is  not  graphically  precise,  it  does 
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provido  tho  user  with  a baso  from  which  ho  can  mako  intu- 
itive interpretations  of  tho  distribution  as  it  relates  to 
his  data. 

A separate  subroutine  is  used  to  read  in  data  to 
SXMFIT.  This  anablos  tho  rosoarchor  to  input  data  from 
cards,  tape,  or  a time-sharing  file  with  a minimum  of 
data  manipulation.  The  data  may  bo  road  in  either  a 
formatted  or  froo-fiold  mode  depending  on  tho  dosiros 
or  noods  of  the  usor. 

The  program  incorporates  tho  flexibility  for  tho 
usor  to:  (1)  specify  tho  parameters  of  the  distribution 
ho  wishos  to  test,  or  (2)  allow  SXMFIT  to  estimate  distri- 
bution parameters  from  tho  data.  Tho  former  option  allows 
the  usor  to  oxamino  his  data  with  respect  to  a specific 
member  of  a family  of  distributions.  When  tho  latter 
option  is  soloctod,  care  must  bo  taken  to  note  that  uso 
of  the  K-S  .tost  may  not  provide  reliable  rasults  when 
population  parameters  aro  estimated  from  sample  data 
(12:86).  Tho  x1  G-O-F  tost  is  alwayo  used  ns  if  tho 
parameters  of  tho  distribution  havo  boon  estimated 
yiolding  a x*  statistic  which  indicates  a loss  of  ono 
dogreo  of  f random  for  oacn  distribution  parameter. 

Rocont  Modifications  to  SXMFIT 

In  January  1977,  tho  RAND  Corporation  completed 
a modification  to  SXMFIT.  Several  changes  were  made  to 
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the  original  program  which  include  tho  incorporation  of 
two  distributions  not  in  earlier  versions:  tho  negative 
and  positive  binomial  distributions.  Other  modifications 
were  mado  to  facilitate  tho  input  of  data,  format  the  out- 
put in  a more  usable  form,  and  correct  some  of  the  compu- 
tational problems  which  had  boon  identified  in  earlier 
usage  of  tho  program,  Tho  result  is  a flexible,  user- 
oriented,  rapid,  and  relatively  easy  to  use  tool  for 
analyzing  tho  distribution  of  largo  or  small  sets  of 
data. 
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